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We derive a generalized form of the Electric Dipole Spin Resonance (EDSR) Hamiltonian in the 
presence of the spin-orbit interaction for single spins in an elliptic quantum dot (QD) subject to 
an arbitrary (in both direction and magnitude) applied magnetic field. We predict a nonlinear 
behavior of the Rabi frequency as a function of the magnetic field for sufficiently large Zeeman 
energies, and present a microscopic expression for the anisotropic electron g-tensor. Similarly, 
an EDSR Hamiltonian is devised for two spins confined in a double quantum dot (DQD), where 
coherent Rabi oscillations between the singlet and triplet states are induced by jittering the inter- 
dot distance at the resonance frequency. Finally, we calculate two-electron-spin relaxation rates 
due to phonon emission, for both in-plane and perpendicular magnetic fields. Our results have 
immediate applications to current EDSR experiments on nanowire QDs, g-factor optimization of 
confined carriers, and spin decay measurements in DQD spin-orbit qubits. 

PACS numbers: 



An electron spin confined in a semiconductor quan- 
tum dot is a promising candidate for spin-based quantum 
information processing 1 2 . To implement reliable quan- 
tum gates and quantum registers, a complete theoretical 
knowledge of all single- and two- spin- qubit decay chan- 
nels is essential in order to mitigate the decoherence and 
relaxation of qubits. Spin-orbit interaction is one of the 
most important mechanisms of spin mixing of different 
orbital states in semiconductor structures. When com- 
bined with phonons or charge fluctuations, this spin mix- 
ture results in a variety of spin decay channels, although 
they have all turned out to be quite weak for confined 
electrons^Ml, leading to impressively long single-spin re- 
laxation time (Ti) (of the order of seconds at 1 Tesla in 
GaAs and minutes in Si:P)P^l The dominant electron 
spin decoherence channel for GaAs QD and Si donors 
is pure dephasing due to hyperfine interaction with the 
environmental nuclear spins^, with the corresponding 
coherence time (T2) having been pushed to few hundreds 
of microseconds in GaAs 8 and a few seconds in Si:F^. 
This significant separation of T\ and T2 time scales means 
that spin-orbit interaction can play a constructive role in 
qubits with a relatively strong spin-orbit interaction. 

Exploiting the orbital part of the electrons and holes 
to manipulate their spins (via the spin-orbit coupling) is 
becoming a common practice in recent experimental and 
theoretical works on extended and confined elect rons in 
heterostructures with strong spin-orbit interact iorPDMUl 
GigaHertz manipulation of confined electrons in semicon- 
ductor QDs is achievable, in the presence of an applied 
magnetic field, accompanied by an ac electric fielcH For 
small Zeeman splitting (compared to the orbital exci- 
tations), the Rabi frequency shows a linear dependence 
on both applied magnetic and electric fields. This spin- 
electric-field coupling was observed in transport experi- 
ments on lateral GaAs QDs and In As nanowires in the 
spin-blockade regime^"^, where the Rabi spin-flip time 
of 10 ns has been reported. Such a fast spin rotation is 



a key step towards realizing single- and two-qubit gates 
for localized spin qubitf^l 

In this paper, we first present the general EDSR Hamil- 
tonian for a single electron spin confined in an elliptic 2D 
QD, in the presence of an arbitrarily large magnetic field. 
The study of elliptic confinement is crucial in nano-wire 
QDs, where the electron o rbital wave function is strongly 
squee zed in two direction J 16 * 17 l In contrast to previous 
work&SMl, we show that the Rabi frequency exhibits a 
pronounced non-linear behavior as a function of the ap- 
plied magnetic field for sufficiently large QDs and Zee- 
man energies. We also calculate the second-order cor- 
rections of the spin-orbit coupling to the EDSR Hamil- 
tonian, and provide an analytical microscopic expression 
for the g-factor renormalization. The ability to engineer 
the g-factor of a confined electron is essential for individ- 
ual addressing of spin-orbit qubits via electric pulses^HH. 
In addition, we propose a scheme to perform coherent ro- 
tation in the singlet-triple subspace of two electron spins 
confined in a double QD. Again, Zeeman energy and the 
spin-orbit interaction play the central role in the electri- 
cal manipulation of the two- spin states, where the corre- 
sponding Rabi frequency is linear, in the leading order, in 
magnetic field and the spin-oribt coupling. The proposed 
two-electron EDSR offers an alternative way to imple- 
ment coherent swap gates in coupled spin-orbit qubits, 
without a magnetic field gradient. Finally, we calculate 
the phonon-induced spin relaxation time of the two-spin 
system and find that the resulting time scales are rela- 
tively long, which justifies the efficiency of our two- spin 
EDSR scheme. 

The effective EDSR Hamiltonian for an elliptic dot. 
Consider a single electron in an elliptic 2D quantum dot, 
with confining frequencies uj x and uo yi subject to an ap- 
plied magnetic field B = £?(sin# cos </>, sin#sin(/>, cos#). 
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TABLE I: The {cti,ai} coefficients of the transformation matrix S in Eqs. ( 5|7 ) as a function of the Zeeman interaction, 
the spin-orbit lengths, the cyclotron frequency uj c = eB z /mc, and the lateral confining frequencies. Note that M is a com- 
mon factor which appears in all coefficients. Similarly, {A? A} are given by the following substitutions: 0i — ai(cj x — »• 
ujy,X- — ► A+), 02 — —ot2{w x — > Wy,\- — > A+), /3i — — ai(A- — > A+), and 02 — oli(w x — ► w y ,\- — > A+); For example, 
02 = hEz(huj c )(hw y ) 2 /mM\+. 
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FIG. 1: (color online). The Rabi frequency versus the mag- 
netic field for an electrically driven single spin in a typical 
InAs nanowire QEM^ uj x — 3 meV, u y — 40 meV, and 
A+ rsj 100 nm, for both perpendicular (ojr±_) and in-plane 
(ujr\\) magnetic fields. The electric field is in-plane and along 
the x axis. The non-linear dependence shows up for B fields 
above 2 T, although this threshold can be substantially re- 
duced for larger dots and/or for materials with stronger spin- 
orbit interaction. 



The Hamiltonian for the electron is 

H = H d + H z + H so , (1) 

H d = ^ + \m^+uy\ (2) 

H z = ^9I^bB ■ a = ^E z n ■ <r, (3) 

H so = 1C ■ a. K = ^{p y /\-,p x /\+,G), (4) 

where m is the effective mass of the electron, p — 
—ihd/dr + (e/c)A(r) the electron kinetic momentum in 
two dimensions, and c the speed of light in vacuum. 
X± = h/m{0 zb a) are the effective spin-orbit lengths 
in terms of the li near Rashba (a) and Dresselhaus (0) 
coupling constant J 18 | 19 l We restrict our consideration to 
quantum dots with strong confinement along one axis 
([001]), such as, e.g., quantum dots defined in a two- 
dimensional electron gas (2DEG). With a uniform B 



field, the vector potential is A(r) = B z (—y/2,x/2 : 0) in 
the symmetric gauge. With the motion along z strongly 
quantized (the 2DEG thickness V <C y/hc/eB^^), the 
in-plane components B x and B y are not included in 
A(r). The magnetic field also induces a Zeeman splitting 
Ez — 9PbB^ with a spin quantization axis n = B/B. 
Note that the linear-in-momentum spin-orbit interaction 
is given in a coordinate system which is rotated by 7r/4 
along the z axis with respect to the crystallographic axes 
of the crystal and it is the sum of both Rashba and Dres- 
selhaus terms. 

In order to diagonalize the Hamiltonian in Eq. ([!]), 
we invoke the Schrieffer- Wolff transformation H — >• 
e~ s He^^Ql, and remove the spin-oribit interaction in the 
leading order by imposing [Hd + Hz,S] = H so . The re- 
sulting S matrix is presented here to all orders in the 
Zeeman interaction, but in the first order in spin-orbit 
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The coefficients {c^, ft, c^, 0i} are given in Table I. Due 
to the spin-orbit interaction, the spin and orbital part 
of the rotated states are strongly coupled to each other, 
which opens up the possibility of spin manipulation via 
driving by an ac electric field V(r,t) = eE • r (within 
the dipole approximation). The effective Hamiltonian, 
Hedsr = H z + (ipo\[S, V(r,t))\^ ) for the spin sector 
of an electron confined in the ground orbital state of 
an elliptic dot, is obtained by averaging over the orbital 
ground state 

H EDS R = H z + \h(t) ■ a, (8) 

h(t) = 2n x {F+g\ (9) 

T = eft (ai£ y , 01^,0), 

Q = eht-nzfaEy.rizaiEx.rixfiiEy -n y aiE x ). 



We note that the effective ac magnetic field h(t) van- 
ishes in the absence of the applied magnetic field, and is 
always perpendicular to B . This is a distinct feature 
of the linear spin-orbit interaction in the leading order. 
Moreover, for an in-plane magnetic field, the vector Q is 
identically zero because uj c = (see Table I). 

In general, Rabi frequency is strongly dependent on 
the directions of the magnetic and the ac electric fields, 
the confinement frequencies, and the orientation of the 
QD with respect to the crystallographic axes. For small 
Zeeman energies (Ez <C hu x ,hu y ), the Rabi frequency 
is linear in Ez 12 . In Fig. [l] we have plotted the Rabi 
frequency ujr = max \h(t)\/h, between the Zeeman sub- 
levels, as a function of the applied B field. A clear non- 
linear behavior arises for magnetic fields above 2 T. To 
obtain these results, we used parameters for a typical 
In As nanowire QD. This nonlinearity is a direct conse- 
quence of the exact treatment of the Zeeman interaction 
in calculating the Schrieffer- Wolff matrix S. 

The above consideration is restricted to a two- 
dimensional QD where the effective spin-orbit inter- 
action is linear in momentum. In reality, all QDs 
have finite thickness V of the order of a few nanome- 
ters. Therefore, the cubic Dresselhaus term ~ 

^7p~ (PyPxPy^x — PxPyPx^y) can also, in principle, induce 
spin transitions 18 , though in the leading order it is in- 
effective in the presence of harmonic confinements and 
in-plane magnetic fields^. 

For spin manipulations in QDs, a co mprehe nsive 
knowledge of the electron g-tensor is cruciaP^^. We 
have calculated the second-order corrections of the lin- 
ear spin-orbit interaction to the effective Hamiltonian, 
Eq. (|8]), which lead to the renormalization of the g- 
factoP^I] For an in-plane magnetic field along the x 
axis, the renormalized g-tensor of an electron in the 
ground state of the harmonic potential reads (results for 
an arbitrary B field direction is presented in appendix 
A) 
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This renormalization of the g-factor can be used to ad- 
dress nearby spins selectively, and has immediate ap- 
p lie at i ons in current EDSR experiments on confined 

carriers^ 14 * 15 *. 

Two- electron spin Hamiltonian. We now consider two 
electrons confined in a DQD separated by 2d, as shown 
in Fig. |2j and the vector d denotes the orientation of 
the DQD with respect to the crystallographic axes. This 
setup has been the cornerstone of recent experimen- 
tal works on spin-selective transport, where the lifetime 
of different spin states a re examined through the spin- 
blockade effect! 13 " 15 l 22 l 23 [ In the absence of the spin-orbit 
interaction, the low energy subspace of the two-electron 
Hilbert space consists of the following four states 




FIG. 2: (color online). Schematic of a singlet-triplet DQD 
spin qubilP^ in the presence of an in-plane magnetic field and 
an ac electric field. The distance between the center of the 
double elliptic dots is 2d and the driving electric field is ap- 
plied along the DQD in order to induce maximum electric 
dipole moment. In the upper right, the corresponding energy 
levels of two-spins has been shown for the regime J <C Ez . 



where S and T±,o refer to the spin singlet and triplet 
states, respectively, and \£± are their corresponding sym- 
metric and anti- symmetric orbital wavefunctions. By in- 
cluding the spin-orbit coupling, the singlet and triplet 
states are coupled provided that an external magnetic 
field is applied. For small magnetic fields, the effective 
Hamiltonian for the two spin confined in a double dot is 
given by 2 -^ 



H ex = ( J/4) <ti • a 2 + fiB • (<ri + cr 2 ) 
+ a • (a i - <r 2 ) + b • (a i x cr 2 ) 



(12) 



where \i = g/i#/2 and J is the singlet-triplet exchange 
splitting. Here we keep only the leading order terms in 
spin-orbit coupling and neglect the second and higher 
order terms. The only linear terms in spin-orbit (and 
Zeeman) interaction are the following vectors 



a = -fiB x Re 
b = -fiB x Im 



(13) 
(14) 



where a — \a\ and b = |6| denote the strength of the 
spin-orbit matrix elements. Note that the Hamiltonian in 
Eq. (12) has the rotated spin-orbit states {U°&i}, where 
U° — exp [-S(°)] 24 as the basis. For an in-plane mag- 
netic field, S and To are coupled via the z component of 
the vector a (there is no coupling, in the leading order, 
between these two states for a perpendicular magnetic 
field). Specifically, a z = (in the Heitler-London ap- 
proximation), provided that B and d are along the x 
axis. Rabi Oscillations between the singlet and the triplet 
states have already been observed in DQDs using the hy- 
perfine field gradient 22 . When spin-orbit interaction is 
included in the consideration, a time-dependent a z can 
also induce coherent oscillations between S and To states, 
providing an alternative way to manipulate the singlet 
and triplet states. For example, a sinusoidal electric field 
pulse can change the inter-dot distance periodically to 
form a breathing DQD, d = do + 5sm(ut). The corre- 
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FIG. 3: The logarithmic plot of the triplet-singlet (To-S) re- 
laxation time, for a DQD along the x axis, as a function of 
the distance between the two QDs. The applied magnetic field 
B = IT ( in this case, the rate has quadratic dependence on 
the field) . The inset shows the corresponding dependence of J 
on 2d, calculated in the Heitler-London approximation. The 
material parameters are chosen for GaAs quantum dots with 



A 4 



pm, 



22 nm, e±4 — 1.4 x 10 9 V/m, 



13.1, 



p = 5.3 x 1(T kg/m^, and s t = 2.5 x 1(T m/s. 



where S <C do is the amplitude of the breathing DQD, 
and S is linearly proportional to the applied ac electric 
field E. For InAs QDs, where the g-factor g « 10 and 
A + ~ 100 nm, a 5 ~ 1 nm at B = 1 T would allow a ir/2 
pulse on the order of 1 ns. 



The exchange coupling between the two spins also os- 
cillates as we periodically change their relative distance. 
Geometrically, this leads to a breathing Bloch sphere 
which beats at the driving frequency. However, if the 
amplitude of modulations is reasonably small, they yield 
small deviations 5 J from the equilibrium value J and 
the coherent Rabi oscillations should still be observable 
in the leading order of the spin-orbit interaction^. 



Singlet-triplet relaxation rates. Due to the spin mix- 
ing via the spin-orbit interaction, phonon emission can 
lead to relaxation (and/or leakage out) of the two-spin 
singlet and triplet states. In III-V semiconductors like 
GaAs or InAs, piezoelectric coupling to acoustic phonons 
is the dominant relaxation channeP^. With an in-plane 
magnetic field (in which regime S-T+ and To-T+ do not 
couple in the leading order) and in the leading order of 
the spin-orbit coupling, the triplet-singlet relaxation rate 
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where ei4 and k are the piezoelectric and the dielec- 
tric constants of the material, respectively, p is the den- 
sity of the sample, s t is the speed of transverse acoustic 
phonons, and I is the size of each quantum dot. This 
relaxation rate is valid when a z <C J. Moreover, a z has 
a geometrical structure which enables us to mitigate this 
relaxation channel by carefully choosing the direction of 
the magnetic field with respect to the crystallographic 
axes and the DQD orientation. In Fig. [3j we have plot- 
ted the relaxation rate in Eq. ( fl5] ), together with the 
exchange splitting (within the Heitler-London approxi- 
mation), as a function of the DQD separation 2d in the 
presence of an applied in-plane magnetic field of 1 T. The 
resulting time scales are relatively long for small Zee- 
man energies (below 1 Tesla). Therefore, we conclude 
that in this regime single-spin and/or hyperfine-induced 
triplet-singlet relaxations dominate the spin relaxation 
processes, while dephasing is probably dominated by hy- 
perfine interaction 9 , charge noise 26 , and phonons^. Note 
that while the Hamiltonian in Eq. (12) is quite generic 
for any inter-dot coupling, the analytical results pre- 
sented here are for a relatively narrow inter-dot distance 
regime, because our analytical calculation is bound from 
the small d by the validity of our Heitler-London model 
to calculate J and a zi and from the large d by the re- 
quirement that singlet and triplet states are close to the 
two-electron eigenstates. 

For a perpendicular magnetic field, the only states 
which are coupled, in the leading order of the spin-orbit 
interaction, are S and T + (See appendix B). Although 
there is no analytical form for this relaxation rate, we 
have found numerical values showing very slow leakage 
rate (0.001 s _1 or slower when a, b <C Ez — J) for mag- 
netic fields of interest between 0.1 - 1 T. 

In conclusion, we have derived the most general form of 
the EDSR Hamiltonian for an elliptic quantum dot, in the 
presence of the (linear) spin-orbit interaction. We find a 
nonlinear 5-field dependence for the Rabi frequency of 
the resulting single-spin EDSR. We also propose electri- 
cal manipulation of two-electron-spin states via inter-dot 
modulation. We find that the phonon-induced triplet- 
singlet relaxation rates are typically small, so that coher- 
ent Rabi oscillations can, in principle, be observed with 
current experimental setups in DQD spin-orbit qubits. 

We thank financial support from NSA/LPS through 
US ARO, and DARPA QuEST through AFOSR. 
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Appendix A: The renormalized g-tensor for confined 
electrons 

In this subsection, we calculate the renormalization of 
the confined electron g-factor in the presence of an ap- 
plied magnetic field. Through the Schrieffer- Wolff trans- 
formation, the effective Hamiltonian (up to second order 
in spin-orbit interaction and to all orders in Zeeman cou- 
pling) is [5,12] 



H eff = H d + H z + AH, 
A# = ±[S,H ao ], 



(Al) 
(A2) 



where the transformation matrix S is introduced in the 
main text of the paper. By calculating the above com- 
mutator, one can obtain the (second-order) corrections to 
the Hamiltonian and the electron spin g-tensor. The re- 
sulting Hamiltonian has a rather complicated form, thus 
we only present two special cases of particular interest: 



1. B in-plane and along the x axis 
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where {/7, V} = UV + VU is the anticommutator of U 
and V. In an in-plane B field, generally all components of 
the g-tensor is renormalized, depending on the electron 
orbital state. However, for an electron in the ground 



state, the quantum mechanical average of the orbital op- 
erators multiplying <j y and a z vanish, and only g xx differs 
from the bulk g-factor. 

2. B perpendicular to the x-y plane 



AH i 
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In this case, only g z 
orbital state. 



is renormalized, for any electronic 



Appendix B: The leakage rate out of the 
singlet-triplet subspace 



regime where a, b <C Ez — J- In a perpendicular B field, 
there is no coupling between S — To, in the leading order 
of the spin-orbit interaction. The only leakage channel is 
S -> T + [see Eq. (12) in the text], 



Here we calculate the two-electron spin relaxation rate 
for a double dot in a perpendicular magnetic field, for the 
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where q = V2 r/l, r l(E z - J)/(y/2Hs t ), I 
IbI\J\ + % Ib = yfjj, and u; c = ^. The function 
/(#) is introduced in the main text of the paper [see 
Eq. (14)], and I n {x) are the modified Bessel functions 
of the n-th kind. Note that the rate in Eq. ( |B1[ ) has a 
geometrical dependence on the orientation of the setup 
with respect to the crystallographic axes. In Fig. [4j we 
have plotted this leakage rate as a function of the mag- 
netic field for a DQD along the x axis. We find that 
the rates are typically very small and much slower than 
the single-spin (and the hyperfine-induced two-spin) re- 
laxation rates. 



FIG. 4: The singlet- triplet (S — » T+) relaxation time, for a 
DQD along the x axis, as a function of applied (perpendicular) 
magnetic field. The material parameters are chosen for GaAs 
quantum dots with A so — 5 /mi, / — 30 nm (at D — 0), 
ei4 = 1-4 x 10 9 V/m, k = 13.1, p = 5.3 x 10 3 kg/m 3 , and 
s t 2.5 x 10 3 m/s. 
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